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Abstract 

Several formulas for crossing functions arising in the continuum limit 
of critical two-dimensional percolation models are studied. These include 
Watts's formula for the horizontal-vertical crossing probability and Cardy's 
new formula for the expected number of crossing clusters. It is shown that 
for lattices where conformal invariance holds, they simplify when the spatial 
domain is taken to be the interior of an equilateral triangle. The two crossing 
functions can be expressed in terms of an equianharmonic elliptic function 
with a triangular rotational symmetry. This suggests that rigorous proofs of 
Watts's formula and Cardy's new formula will be easiest to construct if the 
underlying lattice is triangular. The simplification in a triangular domain of 
Schramm's 'bulk Cardy's formula' is also studied. 
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1 Introduction 



The critical behavior of percolation is not fully understood, either rigorously or 
formally. As the percolation threshold is approached, connected clusters occur 
with high probability on ever larger length scales. It has been conjectured that in 
any dimension, there is a universal scaling limit of isotropic short-range perco- 
lation models, defined over the continuum and independent of the details of the 
model, such as the lattice and percolation type (site or bond) [0]. This continuum 
theory would capture the connectivity of typical configurations of the underlying 
discrete model, at or near criticality. 

Conformal field theory makes predictions for the crossing probabilities of the 
continuum limit of critical two-dimensional percolation, and especially, predicts 
that they are conformally invariant [^. For example, if the continuum theory is 
confined to a spatial domain n plus boundary d^l, and 71, 72 are disjoint pieces 
of dil, the probability of the event that 71, 72 are connected by a percolation cluster 
is predicted to be invariant under transformations that are conformal on (though 
not necessarily on d^l). The probabilities of more complicated crossing events, 
involving more than two pieces of dVl, are also predicted to be invariant. 

In the first applications of conformal field theory to percolation, Q was taken to 

be a rectangle, with aspect ratio r =^ width/height. In this geometry, Cardy 
derived a formula for the crossing function n/j(r), the probability that the two 
vertical sides are horizontally connected. His formula takes on a simpler form 
if the rectangle is conformally mapped onto the upper half plane EI C C, and 
its boundary to dM = M U {00}. The vertical sides are mapped to disjoint line 
segments on the real axis, one of which can be taken without loss of generality to 
be semi-infinite. They are usually taken to be [0, z] and [1, 00], where z E (0, 1), 

with z = 0, 1/2, 1 corresponding to r = 00, 1, 0. We write %^h{z) =^ Uh {r{z)). 
If the underlying discrete model is bond percolation on a square lattice, duality 
suggests Uhil/r) = 1 - Uhir), i.e., ^^,(1 - z) = 1 - ^^iz), so that n;,(l) = 
*P/i(l/2), the probability that two opposite sides of a large square are connected 
by a critical percolation cluster, should equal 1/2. 

On the numerical side, crossing events and their conformal invariance were 
extensively investigated by Langlands and collaborators [1T3|,[T3], and it was veri- 
fied that Cardy's formula is valid for discrete percolation models on a rectangular 
square lattice of size Lx L', with r = L/L', in the limit L, L' 00. They also in- 
vestigated Ilh^(r), the probability that all four sides of the rectangle are connected. 



Watts [ p3| ] derived a formula for the equivalent function ^hv{z) == Hhv {r{z)) 
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from conformal field theory by making additional assumptions, and his formula 
agrees well with the data of Langlands et al. It should be noted that < 
and that by symmetry, Uhvil/r) = nto(r), i.e., ^3^(1 - z) = ^hviz). Also, it is 
clear that Uhv{r) /Uh{r) ^ 1 as r oo, i.e., ^hv{z) /%ih{z) ^ 1 as 2: ^ 0. 

Cardy [|], ^ later derived a formula for the expected number of percolation 
clusters that cross between 71, 72 (the left and right-hand sides of a rectangle, or 
the two corresponding line segments in dW). This may be viewed as a function 

Nh{r), or equivalently a function ''Jlhiz) =^ Nh {r{z)); necessarily, D^/^ > '^h- 

So in all, three crossing formulas for critical percolation have been derived 
from conformal field theory: formulas for ^hviz), and ^h{z), which yield 

formulas for Hhir), Ilhv{r), and Nh{r). Recently, Smirnov p2| ] provided the first 
rigorous proof of any of these, namely Cardy 's formula for ^^(z). He showed that 
critical site percolation on the triangular lattice has a conformally invariant scal- 
ing limit, and that discrete percolation cluster boundaries converge to a stochastic 
Loewner evolution process. These facts followed from his proof that in a general 
domain fi, a conformally transformed Cardy's formula is valid. This includes a 
version due to L. Carleson (unpublished, but see ref. ^ sec. 7.2, and ref. predic- 
tion 7). Suppose HI is mapped conformally onto an equilateral triangle AABC in 
the complex plane, and that the map extends to the boundary, with [—00, 0], [0, 1], 
[1, 00] being mapped to the edges AB, EC, CA. So w e BC and z G (0, 1) will 

correspond. Carleson noticed that ^hiw) == (ziw)) simply equals Bw/BC. 
That is, Cardy's formula is formally equivalent to the statement that that the line 
segment Bw is connected by a critical percolation cluster to the opposite side of 
the triangle, CA, with probability equal to the fraction of the side BC occupied 
by Bw. Equivalently, the real- valued function is the restriction to one side 

of the triangle of a particularly simple analytic function ofw.a linear function. 

The triangular symmetry of Carleson's restatement made possible Smimov's 
proof, which is specific to a triangular lattice. Smirnov notes, "It seems that 27r/3 
rotational symmetry enters in our paper not because of the specific lattice we 
consider, but rather [because it] manifests some symmetry laws characteristic to 
(continuum) percolation." Whether his proof extends to other lattices is unclear. 

In this paper we study whether the predicted formulas for the functions '^hv{z) 
and ^h{z), like Cardy's formula for ^h{z), simplify when the spatial domain Vt is 
taken to be an equilateral triangle, rather than a rectangle or the upper half plane H. 
We show that they do. In particular, we show that the four-way crossing func- 
tion ^h^{w) =^ {z{w)) predicted for the equilateral triangle has its second 
derivative ^'^^^(w) equal to a familiar elliptic function: an equianharmonic Weier- 
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strass p-function, where 'equianharmonic' signifies that the period lattice of the 
p-function is triangular, with a 7r/3 rotational symmetry. (See ref. |l], sec. 18.13.) 
This contrasts with the linear function ^h{w), the second derivative of which is 
zero. Our new representation for ^hviw), i.e., for ^hviz) or Uhvir), immediately 
yields a simple closed-form expression for *P/i.u(l/2) = Ilhv{l), the probability 
that all four sides of a large square are connected by a critical percolation cluster; 
namely, 1/4 + (\/3/47r)(3 log3 - 4 log 2) ^ 0.322. It would be difficuk though 
not impossible to derive this expression directly from Watts's formula. 

We show that Cardy's recent formula for ^hiz), like Watts's formula for 
%ihviz), simplifies in an equilateral triangle. the second derivative of 

=^ ^{z{w)), can also be expressed in terms of the equianharmonic p- 
function. In fact, we derive a curious identity relating the three crossing func- 
tions ^h, ^hv, and or equivalently %^hv, and D^^; namely, that 2^hiz) — 
%ihiz) - ^hviz) must equal (V3/27r) log (1/(1 - z)). Setting z = 1/2 yields that 
D^/i(l/2), i.e., Nh(l), the expected number of critical percolation clusters crossing 
between two opposite sides of a large square, should equal 3 /8+ (\/3/87r) (3 log 3— 
2 log 2) ^ 0.507. In a final study of a crossing formula, we treat a fourth Cardy- 
type formula proved rigorously by Schramm [ [T9| ] for triangular-lattice site perco- 
lation, which is valid 'in the bulk' and does not follow from boundary conformal 
field theory. We show it has a simple restatement in a suitable triangular domain. 

Our successful simplification of the crossing event formulas indicates that an 
equilateral triangular domain is a good 'fit' to the continuum limit of critical per- 
colation. It also suggests that rigorous proofs of Watts's formula and Cardy's new 
formula will be easiest to construct if the underlying lattice is triangular. Our 



restatements contrast with those of Ziff [|28|, [2 90 and Kleban and Zagier [ |11| , |12| ], 
which also involve higher transcendental functions. They focused on the contin- 
uum limit of percolation in a rectangle, and especially on the derivatives Il'h{r) 
and H'fivii^) of the rectangular crossing functions. In our notation, Ziff showed 
that n'^(r) is proportional to [^^'i(0, q = e~'^^')]^^'^ , where 9i{-, q) is the first Jacobi 
theta function. Also, Kleban and Zagier showed that 11'^, 11'^^,, considered jointly, 
have interesting modular transformation properties, and that these properties char- 
acterize n'^^(r). The relation between our results and theirs is not yet clear. 

Section ^ presents each crossing event formula in a standard form. Section |^ 
covers conformal mapping concepts, including the equianharmonic p-function. 
The simplified versions of the crossing formulas that apply in triangular domains 
are derived in Section ^. An appendix reviews some basic mathematical facts. 
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2 Crossing Event Formulas 



The four crossing formulas for continuum percolation on the closed upper half 
plane EI involve hypergeometric functions, both Gauss's 2F1 and Clausen's 3F2. 
They can be stated in a standardized, P-symbol form. (For hypergeometric func- 
tions and P-symbols, see the appendix.) For the first three formulas, dM is divided 
into [—00,0], [0,z], [z,l], and [l,oo]. A connection between [0,^] and [I, 00] 
corresponds to a horizontal crossing on the original rectangle, and one between 
[—00, 0] and [z, 1] to a vertical crossing. The probability of a horizontal connec- 
tion is ^h(z), with ^h{z) the expected number of such connections. All four 

segments are connected with probability Let =^ — ^hv, the 

probability of there being a horizontal connection that is not also a vertical one. 

Formula 2.1 (Cardy [0]) The function ^hiz) equals 
3r(2/3) , 
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Formula 2.2 (Watts [^) The function ^^hv{z) equals 
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Formula 2.3 (Cardy The function ^h{z) equals 

1, 1, 4/3 



1 

2 ~ 4^ 



loEcfl - z) 



'l-z) 3F2 



2, 5/3 



1 - z 



(3) 



Proof. This version of the formula for ^h{z) is not well known. The version 
deduced by Cardy from boundary conformal field theory was 



1 V3 

2 ~ 4^ 



log(l-2;)+2 5^ 



(1/3)^, (1 - zy 

\ (2/3)m, 



m 



(4) 



Formula ^3] follows from the series representation (|A1[), if the 1/m factor in the 
summand is written as (l)m_i(l)m_i/(2)m_i(m — 1)!. ■ 
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Corollary 2.3.1 2mh{z) - ^h{z) - ^/,.(^) = (v^/27r) log (1/(1 - z)). 

Proof. This follows by combining Formulas p]-|2.3|, with the aid of the symmetry 
relations ^^{1 - z) = I - ^h{z) and ^^vil - z) = ^h,{z). ■ 

Corollary |2.3.1| ties to and ^P/j^, in quite a strong way. To explain how, 
we must sketch the heuristic origins of Formula in boundary conformal field 
theory. Watts was led to Formula by considering ODEs satisfied by correlation 
functions of boundary operators. His candidate for an ODE satisfied by "^hv was 
the fifth-order Fuchsian equation 
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which has P- symbol 
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but is not of hypergeometric type. Due to a factorization of the differential op- 
erator on the left-hand side of this equation [[TT],p^, its solution space properly 
contains the solution space of the third-order equation 



0, 



(7) 



which has P- symbol 
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and is of hypergeometric type. Furthermore, the solution space of properly 
contains the solution space of the second-order equation with P- symbol 
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including the function of Formula P7T] . Watts noticed that the solution space 
of includes a one-dimensional subspace of functions equal to zero at 2; = 
and invariant under z t-^ 1 — z, which are criteria for the function %ihv For this 
reason, he expected ^hv^ and %ihv = — %^hv too, to be solutions of (^, as well 
as of This insight led to Formula |Z2|, which incorporates the P-symbol 

Kleban and Zagier [ [T^ noticed that the five-dimensional solution space of @ 
is spanned by the solution space of (^, to which ^^.^ and belong, and the 
functions log z and log(l — z). But by Corollary |2.3.1| , the function is a linear 
combination of ^hv, ^h, and log(l — z), implying the following mysterious fact. 

Corollary 2.3.2 The function ^hiz), like the functions %ihviz) and ^hiz), is a 
solution of Watts's fifth-order differential equation, Eq. (|^). 

For the fourth crossing formula, dM is divided into [—00, z] and [z, 00], with 
2 G M unrestricted. A special boundary condition is imposed: on the underlying 
discrete lattice, percolation along the line segment [—00, z] is allowed by fiat. Let 
a distinguished point in EI be chosen; without loss of generality, let it be i = 
Then the function ^smriz) is defined to be the probability that i is surrounded by 
the percolation hull of [—00, z], i.e., the outermost boundary of the percolation 
cluster in EI containing ("growing from") [— cxd,^]. One expects ^smri—z) = 
1 — ^surr{z), i.e., that ^surriz) — 1/2 is odd in z. 



Formula 2.4 (Schramm [|T9p) The function %isurr{z) equals 1/2 plus 
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In the P-symbol expression on the right, z < and ^ > correspond to different 
branches, which are negatives of each other. 



3 Conformal Maps and Fuchsian ODEs 

The upper half plane EI and any triangle AABC without boundary are homeo- 
morphic as complex manifolds. In fact, by the Poincare-Koebe uniformization 

theorem, the complex plane C, the Riemann sphere CP^ =^ C U {oo}, and EI are 
the only simply connected one-dimensional complex manifolds, up to conformal 
equivalence [j^]. To transfer the Fuchsian ODEs satisfied by the crossing func- 
tions from EI to AABC, an explicit expression for the conformal map w = s{z). 
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i.e., s : H — > A ABC, is useful, s is a Schwarz triangle function, defined by 
the Schwarz-Christoffel formula [[18]]. It extends to a map from the closure of H 
in CP^, i.e., HUMU {oo}, to the triangle with boundary. The conditions s(0) = B, 
s(l) = C, s(oo) = A uniquely determine s. 

The inverse Schwarz function S : AABC H often has a deeper signifi- 
cance than s does, as the equilateral triangle case illustrates. Equilateral triangles 
tile the plane, and the inverse function extends to a function S : C ^ CP^ that 
maps alternating triangles in a checkerboard fashion to the upper and lower half 
planes. In modem treatments, this extended map is classed as one of a handful of 
'universal' branched covers of CP^ by C, CP\ or M (see ref. sec. 6.4). 

Since the triangular lattice is doubly periodic, one expects that S is an elliptic 
function, and can equally well be viewed as a function on a complex elliptic curve 

C/£, where C =^ 2ujZ + 2uj'Z is an appropriate lattice of periods (the factors of 2 
are traditional). This is correct, as we briefly sketch; for details, see Abramowitz 
and Stegun [[Ip, sec. 18.13, and Sansone and Gerretsen [[r8|], sees. 14.2-14.3. The 
canonical elliptic function on C is the Weierstrass function p{w; g2, g^), defined 

as the solution of {p'Y = f{p) =^ 4:p^ — g2p — g3 with a unit-strength double pole 
at 1(7 = 0. The parameters g2,g3 G C, both of which cannot be zero, are related in 
a nontrivial way to the fundamental half-periods uj,lu' of p. The equianharmonic 
case, which has special symmetries, is the case when (72 = 0. Due to the homo- 
geneity relation p{w, g2, gs) = t^p{tw; t~^g2, t'^g^), in the equianharmonic case 
all nonzero g^ are equivalent, so henceforth g^ E M \ {0}, in particular g-^ = 1, 
will be taken. In this case the fundamental half-periods u, u' can be chosen to be 
a complex conjugate pair with 'Situ > 0, '^cu < 0. If the basic real half -period 
uj + Lu' > is denoted uj2, then cu, cu' will equal =1= -^^)co'2. So the period lat- 
tice C will be a triangular lattice, with a 7r/3 rotational symmetry about the origin. 
Explicitly, CU2 = r(l/3)V47r ^ 1.530. 

The elliptic curve C/£ is homeomorphic to a torus and can be viewed as the 
parallelogram with vertices 0, 2uj, 2uj2, 2uj', equipped with periodic boundary con- 
ditions. The p-function maps this parallelogram doubly onto CP^. Also, p' maps 
it triply onto CP^. It turns out that in the equianharmonic case, the torus, i.e., 
this period parallelogram, can be subdivided into six equilateral triangles, each 
mapped by p' with unit multiplicity onto the left or right half plane. (See ref. [1|, 
fig. 18.11, which is unfortunately not quite to scale.) Due to this, the equilateral 
inverse Schwarz function S can be chosen to be essentially p'. The map 

z = S{w) = 1/2 + p'{w)/2i (10) 
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del 

will take the equilateral triangle AABC = AOWqWo in w-space to the upper 
half plane in 2;-space, where Wq =^ (1 + ^^^)uj2- It will take the boundary of 
AABC to M U {c>o}, and the vertices w = A, B,C respectively to ^ = oo, 0, 1. 
As a map from C to CP^, it will take alternating triangles to the upper and lower 
half planes. Each of these equilateral triangles, which tile the plane, has one vertex 
in each of the congruence classes A + C,B + C,C + C, i.e., in each of the sets 
S~\oo) = C,S-\0),S-\1). These classes will be denoted [A], [B], [C]. 

The ODEs on CP^ D H that are satisfied by the functions of Formulas p7T|-p3| 
can be pulled back to ODEs on C via the extended map S : C ^ CF\ The 
important thing to note when performing the puUback is that 5 is a branched 
cover of CP^ by C, the critical points of which are the points in [A\, [B], and [C]. 
At any wq G C, necessarily S{w) ~ S{wq) + const x (ty — wqY to leading order, 
where p is the multiplicity with which wq is mapped to S{wq). (If S{wq) = oo, 
the right-hand side must be replaced by const x (w — wq)^'^.) Each critical point 
has p = 3. 

The pulled-back ODEs are tightly constrained by the following lemma, which 
specifies how P-symbols are pulled back. It is proved by considering the local 
(Frobenius) solutions at each point. 

Lemma 3.1 Let R : M ^ CP^ be a holomorphic map of one-dimensional com- 
plex manifolds. Consider Lu = 0, an nth-order Fuchsian ODE on CP\ It can be 
pulled back via i? to a Fuchsian ODE Ltt = on M, in the sense that iiu = u{z) 

satisfies Lu = then u = u{w) =^ u{R{w)) will satisfy Lu = 0. The n char- 
acteristic exponents of L at each G M will equal those of L at zq = R(wo), 
multiplied by the multiplicity with which wq is mapped to zq. 

As an illustration of the use of this lemma in pulling back ODEs, we give a 
new proof of Whipple's quadratic transformation formula for 3F2. (The formula 
originally appeared in ref. |2^, with a combinatorial proof; a simpler combinatorial 
proof is due to Bailey [^. For a useful discussion placing the formula in context, 
see Askey [Q], but note the misprint in eq. 2.19.) This new proof resembles Rie- 
mann's concise P-symbol proof of Kummer's quadratic transformation formulas 
for 2F1, which is summarized in ref. ^ sec. 3.9. It relies on the expression ( |A4l ) 
for the P-symbol associated to q+iFg, which is apparently not well known. 

Proposition 3.2 Let a,b,c G C with neither a — b + I nor a — c + 1 equal to 



9 



a nonpositive integer. Then 

= (1-W)-"3F2 

holds in a neighborhood of w = 0. 
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Remark. The two sides are defined for all w such that \ w\ < 1, resp. for all w such 
that \—4:w/{l — wY\ < 1. So by analytic continuation, equality holds at all points 
within the loop of the curve |4w| = 1 1 — wf surrounding the origin. 

Proof. The functions of w on the two sides satisfy third-order Fuchsian ODEs 
that follow from the q = 2 case of (|A2[), the ODE satisfied by q+iFg. The left 
and right-hand functions are determined by the additional condition that they be 
analytic at ui = and equal unity there. It will therefore suffice to prove that the 
Fuchsian ODEs corresponding to the two sides are the same up to normalization, 
i.e., have the same solution spaces. A necessary condition for this is that their 
P-symbols be the same, i.e., by the representation (|A^, that 
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where R : CF^ ^ CP^ is defined by R{w) = -Aw /{I — w)^, or equivalently 
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That is, R must pull back the right-hand P-symbol in (|T^ to the left-hand one. 

The map i— > z =^ R{w) takes w = 0, 1, cxo to z = 0, oo, respectively, 
and also w = —1 to z = 1. Its critical points are w = ±1, each of which has 
double multiplicity. The exponents in the columns of agree precisely with 



what Lemma pTT] states: the exponents of w = and w = oo are the same as those 
of 2 = 0, and those of w = 1 are twice those of 2; = 00. One might think the 
left-hand P-symbol would have a fourth critical point, atw = —1, with exponents 
twice those of 2; = 1, i.e., 0,1,2. But as noted in the appendix, those exponents 
are the signature of an ordinary point; so no fourth column is present. 

The preceding argument shows why the 3F2 parameters of the proposition take 
the values they do, but it does not quite prove the proposition. As reviewed in the 
appendix, any Fuchsian ODE on CP^ that has three singular points and is of third 
order (i.e., q = 2) has 3g + 2 = 8 independent exponent parameters, which are 
displayed in its P-symbol, and (2) = 1 accessory parameter, which is not. For 
the two ODEs to be the same up to normalization, they must have the same P- 
symbol, and also the same accessory parameter. The latter is most readily verified 
by changing variables from z = R{w) to w in the right-hand ODE. An explicit 
computation, omitted here, shows that the resulting pulled-back ODE is indeed 
the q = 2 case of the ODE (|A2|), with the parameters of the left-hand side. ■ 

What will be used in Section ^ is the following variant of Whipple's quadratic 
transformation. It seems not to have appeared in the literature. 

Proposition 3.3 Let a,b,c G C with {a + b + c)/2 not equal to a nonpositive 
integer. Then 



f a,b, c 

3^2 



w 



2, (a + 6 + c)/2 

/ (a + l)/2, (6 + l)/2, (c+l)/2 
= ^^"^^^^^1 2, (a + 6 + c)/2 



4w(l - w) j 



holds in a neighborhood of w = 0, provided one of a, b, c equals unity. 

Remark. The two sides are defined for all w such that \w\ < 1, resp. for all w 
such that |4w(l — w)\ < 1. So by analytic continuation, equality holds at all 
points that are both within the circle \w\ = 1 and within the loop of the curve 
\4:w{l — w) \ = 1 surrounding the origin. 

Proof. This closely follows that of Proposition |^; the details are left to the 
reader. The only new feature is that equality between the accessory parameters 
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of the ODEs satisfied by tlie two sides leads to an additional condition on the 3F2 
parameters, beside the exponent conditions of Lemma p7T] . Changing variables 

from z = R{w) =^ Aw{l — w) to w in the right-hand ODE pulls it back to the 
left-hand ODE, plus an extraneous term proportional to (a — — l)(c — 1). 
Provided one of a, 6, c equals unity, this undesired term is absent. ■ 



4 Transformed and Restated Formulas 

In this section we show how the crossing event formulas. Formulas p7Tl-|2.4|, sim- 
plify in appropriately chosen triangular domains. The restated formulas appear 
in Propositions ^"1]-^^^. For all but Schramm's formula, the appropriate triangle 
is equilateral. The restatements of Cardy's formula and Schramm's formula do 
not involve special functions. The restatements of Watts's formula and Cardy's 
new formula do involve elliptic functions, but elliptic functions are significantly 
more familiar than Clausen's 3F2. The restatement of Cardy's formula is of course 
identical to Carleson's, but the other three are new. 

The Fuchsian ODEs on CP^ D H of Formulas O-O are pulled back to 



ODEs on C D AABC via the inverse Schwarz function 5 : C ^ CP\ In the nor- 
malization of the last section, z = S{w) equals l/2 + p'{w)/2i, with ptheequian- 
harmonic p-function, satisfying {p'Y = — 1. The equilateral triangle AABC 

is AOW^Wo, i.e., AO, pe-^'^/^ pe''^/^ with side length p = \Wo\ = 2uj2/\^. 

As noted, UJ2 *= r(l/3)^/47r ~ 1.530 is the basic real half-period of p. 

Since S maps AB, BC, CA to [— oo,0], [0, 1], [l,oo] respectively, the line 
segment BC = WqWq is of primary interest. Its midpoint is uj2 = (Wq + Wo)/2, 
which is mapped to 1/2. As a necessary preliminary, the behavior of the first few 
antiderivatives of p' along BC will now be described. (See the tables in ref. [1], 
sec. 18.13, where Wq is denoted '-Zq'.) Relative to the midpoint, p' is an odd 
function: it equals —i at Wq and i at Wq. Its antiderivative p is even: it equals 
zero at Wq and Wq, and 4^^/^ at the midpoint. The negative antiderivative of p 
is the so-called Weierstrass zeta function, plus an arbitrary constant. The shifted 
negative antiderivative ( — -K/2y/3uj2 is odd: it equals m/Qu2 at Wq and —m/Qu2 
at Wq. The antiderivative of C, equals logo" plus an arbitrary constant, where a is 
the Weierstrass sigma function, which equals g'^/Sx/Ig-jTr/e ^-K/zV^^nr/Q 
at Wq, and e''/4^2^/^3~^/^ at the midpoint uj2. It is easily checked that 
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which is a double antiderivative of p, is an even function relative to the midpoint: 
it equals zero at Wo and Wq, and (1/12) (7r/v^+ 3 log 3 - 4 log 2) at the midpoint. 
Log signifies the principal branch of the logarithm function. 

Proposition 4.1 (Cardy's formula, transformed; cf. Carleson) If conformalin- 
variance holds, Formula ^]T] corresponds on the equilateral triangle /\ABC plus 
boundary to the following, ^/^(w), the probability that the boundary segments 
Bw and CAare connected by a percolation cluster, is the restriction to BC of an 
analytic function that is linear. Explicitly, = {w — B) / (C — B). 

Proof. By Lemma |?7T|, the P-symbol of Formula ^7T] is pulled back via z = S{w) 
to 



[A] [B] [C] 



w 



(14) 



p<; 

1 □ 1 

where [A], [C] are the classes of points on the w-plane that are mapped by S 
to z = cxD, 0, 1 respectively. This is because these points are the critical points 
of 5, and each has triple multiplicity. But a singular point with exponents 0, 1 
is effectively an ordinary point. So the pulled-back ODE on C (in particular, 
on /\ABC) has no singular points and should be effectively {d? / dw'^)'^h{w) = 0, 
as can be verified by an explicit computation. ■ 

Proposition 4.2 (Watts 's formula, transformed) If conformal invariance holds. 
Formula corresponds on the equilateral triangle A ABC plus boundary to the 
following, ^hviw), the probability that all four boundary segments AB, Bw, 
wC, and CA are connected by a percolation cluster, is the restriction to BC of 

an analytic function with the property that the difference ^hjj{w) =^ '^hi'w) — 
^hv{w) is proportional to {w — B)^ as w ^ B, to leading order. Explicitly, 

^ / X 3v3 ^ , . 3 w 1 

"^hviw) = hoga{w) + , (15) 

where o is the equianharmonic Weierstrass sigma function. 

Proof. By Lemma PTT] , the P-symbol of Formula P?T| which partially specifies the 
ODE satisfied by ^^^{z) and '^h{z), is pulled back via z = S{w) to 



P 



[A] 


[B] 


[C] 


w 













1 


1 


1 









3 





(16) 
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The condition Uh^^r) /Ilh{r) — > 1 as r — > oo, i.e., %^hy{z)/^h{z) — 1 as 2; — > 0, 
implies ^hviu^)/%^hiw) ^ 1 as w ^ B. So %ihviu^) is linear inwasw^B, 
to leading order. By the P-symbol ([T^), the first nonzero correction must be cubic. 

By changing variables from z = S{w) to w in the third-order ODE satis- 
fied by ^hviz) and %ihiz), one obtains the striking pulled-back ODE 

^{[Piw)]-'VIH]=0 (17) 

on C. So must be proportional to a double antiderivative of p. The con- 
dition Ilhv{l/r) = Ilhv{r), i.e., ^^(1 — 2) = ^hviz), implies that on the line 
segment BC, must be even around the midpoint w = uj2- Moreover, the 
condition ^hviw) /^hiw) ^ 1 as w ^ B implies = at the endpoints 
w = Wq, Wq. Any even double antiderivative ofjo equalling zero atw = Wo, Wq 
must be a constant times the function ([T3[). For ^hv, the constant is set by 

%,{B) = %{B) = 1/iC -B) = l/{Wo - Wo) = v^/2zcu2, (18) 

together with the fact that \Loga\' {B) = (^{Wq), the value of which is given 
above. By examination, the constant should be sVs/n; which yields (|T5|). ■ 

Numerical Remark. A power series expansion of p{w) about w = Wq that is 
accurate to O {{w — Wq)^^) is given in ref. [1], eq. 18.13.41. The corresponding 
expansion about w = Wq, i.e., about B, is obtained by complex conjugation. By 
twice anti-differentiating this, an expansion of ^hviw) about w = B accurate to 
O {{w — By^) can be obtained. 

Corollary 4.2.1 Uhvi^), the probability that all four sides of a large square are 
connected by a critical percolation cluster, equals 

1 /4 + (Vs/An) (3 log 3 - 4 log 2) ^ 0.322. (19) 

Proof. Uhvir = 1) = ^hviz = 1/2) = %ihviw = 0^2) by conformal invariance. 
This quantity can be computed from (jT5|), using the closed-form expression for 
a{uj2) given at the beginning of this section. ■ 

Alternative Proof. The expression ( [T^ ) for Ilhy{l) can be derived directly from 
Watts 's formula, though the derivation is intricate; the following explains how. 
n;,,(l) equals ^^,.(1/2), i.e., Wf^) " ^h^;(l/2). By Formula ^ 

nai) = l/2-^3F.( Y's/'f |l/2), (20) 
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since ^/i(l/2) = 1/2. Summing the 3F2 series requires care, since in general, it is 
harder to evaluate 3F2 than 2F1. For example, though Gauss's formula 



2Fl 



a, b 



T{c)T{c-a-b) 
T{c-a)T{c-b)' 



b)>0, 



(21) 



evaluates any 2-^1 at unit argument, no general formula for 3^2(0, 6, c; d, e; 1) 



in terms of gamma functions exists p^ , |27| ]. However, certain special 3-^2'^ can be 
evaluated at unit argument in closed form. At other argument values, the situation 
is unresolved. Many 'strange' evaluations of 3F2 and 2F1 at rational points other 
than unity are known [|T0|], but most can apparently be reduced to evaluations at 
unity via appropriate transformations of the independent variable. 

To move the 3F2 evaluation point in ( ^0| ) from 1/2 to unity, the new quadratic 
transformation formula of Proposition ^31 can be used. It yields 



n;..(i) 



(22) 



(The point 1/2 is on the boundary of the region to which Proposition ^]3| applies; 
but by the 3fJ(^ /3i — ^ ttj) > convergence criterion mentioned in the appendix, 
the equality of the proposition extends to the boundary.) Fortunately, the 3-^2(1) 
in ( P^ ) can be evaluated in closed form. An ingenious application of L'Hopital's 
rule to Gauss's formula shows that 



1, 1, a 
" 2, c 



3-^2 



(c-l)-V(c-a)], a^l,^{c-a)>0 



(see Luke [|T5]], sec. 5.2.4). Here ip '= F'/F is the digamma function. So 

n,.(l) = l/2-^ [^(2/3)-^(l/2)]. 



(23) 



(24) 



The values ^(2/3), ^(1/2) are -7 + ^l/2^/3 - (3/2) log3 and -7 - 2 log 2 re- 
spectively, where 7 is Euler's constant. (See ref. [T^, vol. 2, app. II.3.) Substitution 
yields the expression ([T^ for 11/1^,(1). ■ 

Proposition 4.3 (Cardy's new formula, transformed) If conformal invariance holds. 
Formula c orresponds on the equilateral triangle A ABC plus boundary to the 
following, ^hiui), the expected number of percolation clusters connecting the 
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boundary segments Bw and CA, is the restriction to BC of an analytic function. 
Explicitly, equals 



V3_( 



6 Loga{w) + Log 



1 P'jw) 

2 2i 



(3 — VSi) w ^/Si 
4 UJ2 A 



Proof. This follows from Corollary |2.3.1| by replacing ^hv{z), z 

by ^/i(u'), *Pto(u'), S{w), respectively, and substituting the expressions 

for ^h^{w) provided by Propositions ^J] and 



Corollary 4.3.1 iV/j(l), the expected number of critical percolation clusters cross- 
ing between opposite sides of a large square, equals 

3/8 + (Ts/Stt) (3 log 3 - 2 log 2) ^ 0.507. (25) 

Proof. Nh{r = 1) = ^hiz = 1/2) = = UJ2) by conformal invariance. 

This quantity can be computed from the formula for 9T/j(w), using the known 
value of cr(ci;2) and the fact that p'{uj2) = 0. More simply, it follows from Corol- 
lary |2.3.1| by substituting the expression for Hhvir = 1) = %ihviz = 1/2) pro- 
vided by Corollary [4.2. 1| . ■ 



Finally we come to Schramm's formula. Formula It differs from Formulas 
P?T|-P3| in that its restatement employs a triangular domain that is not equilateral. 
Let AA'B'C C C be an isosceles triangle with interior angles 27r/3, 7r/6, vr/G 
at A', B', C, respectively. (For concreteness, take the vertices A', B', C equal 
to 1 + iv/3/3, 0, 2, respectively.) Special boundary conditions are imposed: the 
edge B'C is divided into B'w and wC, and on the underlying discrete lattice, 
percolation along B'w is allowed by fiat. Also, for percolation purposes the edges 
A'B' and A'C are identified, so that in effect, the boundary of AA'B'C com- 
prises only the edge B'C, and the vertex A' is in its interior. 

Proposition 4.4 (Schramm's formula, transformed) If conformal invariance holds. 
Formula corresponds on the triangle AA'B'C plus boundary, with edges 
A'B', A'C identified, to the following. Let ^surri^) denote the probability that 
the vertex A' is surrounded by the percolation hull of the boundary segment B'w, 
i.e., the outermost boundary of the percolation cluster growing from B'w. Then 
^surriw) is the restriction to B'C of an analytic function that is linear. Explicitly, 

^surr{w) = {W - B')/{C -B'). 
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Proof. The first thing to observe is that up to trivial changes of the independent 
and dependent variables, the function '^surr{z) of Formula |Z4l is identical to the 
function of Cardy's Formula P?T], or more accurately to its analytic contin- 

uation. This is the source of the linear behavior on AA'B'C To see the close 



relation between the two functions, use Lemma |3J| to pull back the P-symbol of 
Formula via the quadratic map z i-^ —z"^ on CP^. The result of this procedure 
is that ^surr{z) equals 1/2 plus a function in the solution space specified by 



P 






1 


oo 


-z^ 













1/2 


1/3 


1/6 





p 



—i 


i 


oo 


z 













1/3 


1/3 


1/3 





(26) 



since the quadratic map has 0, oo as its critical points (of multiplicity 2), and takes 
0, ±i, oo to 0, 1, oo. There is no fourth column associated to z = in the right- 
hand P-symbol, since in the pulled-back ODE, z = has exponents 0, 1 and is 
effectively an ordinary point. The close connection between this P-symbol and the 
P-symbol of Formula [2. 1| is obvious. A careful computation, omitted here, yields 

^^surriz) = 1/2 + const X [<P/,(l/2 + iz/2) - 1/2] . (27) 



But for the purpose of proving the proposition, ( |2q ) will suffice. 

Via a conformal map R similar to the map S used in the proofs of Propositions 
pTH-|4-3|, the right-hand P-symbol in ( p6| ) can be pulled back to a trivial P-symbol. 



However, it will turn out that R maps a triangle AA'B'C of the above form not 
onto H, but rather onto the slit half plane HI\ [i, +ooi). The edges A'B', A'C will 
be mapped to opposite sides of the slit, and will therefore need to be identified for 
percolation purposes. A' will be mapped to i, so the statements of the proposition 
and Formula |2!4l will correspond. The map R and AA'B'C are chosen as follows. 
The function S{w) = 1/2 + ^'{w)/2i maps AABC = AO, Wq, Wq onto H, 

audits vertices to oo, 0, 1 G = RU{oo}. So R{w) = p'{w) maps AO, W^, Wq 
onto the right half plane, and its vertices to ooi, G Mi U {ooi}. By reflecting 
through the line passing through Wq and Wq (and their midpoint, the real half- 
period UJ2), it follows that as well, R maps the triangle A2c<j2, Wq, Wq onto the left 
half plane. Therefore R maps the parallelogram without boundary 0, Wq, 2ij2, Wq 
comprising these two equilateral triangles and the line segment lioW^o (their com- 
mon boundary) onto the doubly slit plane C \ [i, +ooi) \ {—ooi, —i]. In particular, 
it maps the upper half of this parallelogram, the isosceles triangle AWq, 0, 2uj2, 
onto the slit upper half plane. This triangle has interior angles 27r/3, 7r/6, 7r/6. Its 
vertices are mapped to i, +ooi, +ooi, respectively. 
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One accordingly chooses AA'B'C = AVro,0,2cu2 oc A(l + ^^3/3), 0,2, 
where the constant of proportionality equals 102- The puUback proceeds as in the 
proof of Proposition P~T] . The right-hand P-symbol of is pulled back via 

z = R{w) to 



[A] 


[B] 


[C] 


w 













1 


1 


1 





(28) 



So the pulled-back ODE on AA'B'C (or more generally, on C) has no singular 
points and should be effectively (d'^ / dw'^)^surriw) = 0, as can be verified by an 
explicit computation. Therefore ^smr niust be linear. Since ^smriB') = and 
^surriC) = 1, the proposition follows. ■ 



A Hypergeometric Functions, P-Symbols 

The following are facts about the generalized hypergeometric function q+iFg and 
its ODE [Pj]. Let the rising factorial a{a + 1) ■ ■ ■ (a + /c — 1) be denoted (a)^; 
by convention, (a)o is interpreted as unity. Then for any q > 1, the function 
. . . , aq+i] Pi,..., (3q, z) is defined by 



q+lFq 



A,...,/?, 



I? 



Provided no denominator parameter jSi is a nonpositive integer, the series coeffi- 
cients are finite and the series converges absolutely on the open unit disk 1^1 < 1. 

Provided 3fJ(^ A — X] > 0' it converges on\z\ = 1, as well. Let ?9 =^ z d/dz. 
Then on the disk, q+iFg satisfies the order- (g + 1) ODE 

^9(^9 + /3i - 1) ■ ■ ■ + /3, - 1) - z{^ + ai) ■ ■ ■ + F = 0. (A2) 

It is the only solution analytic at 2; = and equalling unity there. 

The natural domain of definition of ( |A2| ) is the Riemann sphere CP^ =^ CUoo. 
By examination, this ODE has 2; = 0, 1, 00 as its only singular points, and is Fuch- 
sian: each singular point is regular. q+iFg can be continued to a meromorphic 
function on CP^ \ {0, 1, 00}, which is generally multivalued. In fact, the solution 



space of ( |A2D is a (g + 1) -dimensional space of multivalued meromorphic func- 
tions. To avoid multivaluedness, CP^ is cut along [1, 00]. By definition, q+iFg is 
the continuation of the series from the disk to CP^ \ [1, 00]. 
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The solution space of any order-n Fuchsian ODE on CP^ with three singular 
points is determined to a large extent by their locations and the n characteristic 
exponents associated to each []T^]. Let zi, Z2, Z3 G CP^ denote the singular points, 
and p^^\ . . . p'n G C the exponents of z = Zi. If Zi ^ 00, this generally means 

(i) 

that for each j e {1, . . . , n}, the equation has a solution asymptotic to {z — ziY^ 

as z Zi. (If Zi = 00, then {z — zjY^ must be interpreted as z . Also, if 
the difference between any pair of exponents of a singular point is an integer, the 
solution corresponding to the smaller one may include a logarithmic factor.) This 
definition of characteristic exponents extends immediately to ordinary points. The 
n exponents of any finite ordinary point z ^ zi, Z2, Z3 are 0,l,...,n — 1. 
The Riemann P-symbol for such an ODE, or for its solution space, is 





f ^1 


Z2 




\ 

Z 


p< 




(2) 
Pi 


(3) 
Pi 


> 




.pi" 


fJn 




> 



(A3) 



where the order of exponents in each column is not significant. This tableau fa- 
cilitates symbolic manipulation. For example, multiplying the general solution by 
(z — zqY will add c to the exponents of ^ = 2:0 and — c to the exponents of 2 = 00. 

It is readily verified that the ODE ( |A2D has exponents 0, 1 — Pi, . . . ,1 — Pg 
at 2; = 0, exponents 0, 1, 2, . . . , g — 1, s at 2; = 1, and exponents cti, . . . , a^+i 

at 2; = 00, where s =^ ^ /?j — ^ Oj. (This seems not to be well known; it is only 
partially explained in ref. ^.) In an ad hoc notation, we write 



Oil, ■ ■ ■ , Oiq+l 

q 



rp I ^ij ■ ■ ■ J 





f 







1 


00 


z 














ai 




p < 




1 - 


-Pi 


1 


a2 












> 






1 - 


-Pq~l 


q-1 


Oiq 








1 - 


-Pq 


s 


aq+1 


> 



(A4) 



with the box indicating that q+iFg belongs to the zero exponent at 2 = 0. The sum 
of the 3(g + 1) exponents equals (''2^), though this property is not specific to the 
hypergeometric ODE: it holds for any order- (g + 1) Fuchsian ODE on CP^ with 
three singular points. So there are only 3g + 2 independent exponent parameters. 
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Any order-(g + 1) Fuchsian ODE on CP^ with three singular points, or more 
accurately its solution space, is characterized by the 3g + 2 independent exponent 
parameters and (2) additional 'accessory' parameters, which together with the 
exponent parameters determine the global monodromy. (See Poole [|16p, sec. 20, 
for a normal form for the ODE from which the parameters may be extracted, with 
some effort.) The second-order (i.e., q = 1) case is special in that there are no 
accessory parameters, and the solution space of the ODE is uniquely determined 
by its P-symbol. This is not the case when q > 2, i.e., when the ODE is of third 
or higher order. 

If the q + 1 exponents of one of the three singular points are 0, 1, . . . , g — 1, s 
for some s, up to an overall additive constant, we say the ODE is of hypergeomet- 
ric type. Provided its accessory parameters take suitable values, the solutions 
of any ODE of hypergeometric type can be expressed in terms of hypergeometric 
functions, since it can be transformed to the hypergeometric ODE by redefining its 
independent and dependent variables so as to move its singular points to 0, 1, 00, 
and remove the additive constant. 
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